Electron operator at the edge of the 1/3 fractional quantum Hall liquid 
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This study builds upon the work of Palacios and MacDonald (Phys. Rev. Lett. 76, 118 (1996)), 
wherein they identify the bosonic excitations of Wen's approach for the edge of the 1/3 fractional 
quantum Hall state with certain operators introduced by Stone. Using a quantum Monte Carlo 
method, we extend to larger systems containing up to 40 electrons and obtain more accurate ther- 
modynamic limits for various matrix elements for a short range interaction. The results are in 
agreement with those of Palacios and MacDonald for small systems, but offer further insight into 
the detailed approach to the thermodynamic limit. For the short range interaction, the results are 
consistent with the chiral Luttinger liquid predictions. We also study excitations using the Coulomb 
ground state for up to nine electrons to ascertain the effect of interactions on the results; in this 
case our tests of the chiral Luttinger liquid approach are inconclusive. 

PACS numbers: 



I. INTRODUCTION 

The low-energy excitations of an ordinary Landau 
Fermi Uquid resemble electrons, and are perturbatively 
accessible from the free system. In recent years, there 
has been much interest in systems that do not conform 
to the Landau Fermi liquid paradigm. An example of 
a non-Landau Fermi liquid (NLFL) is the system of in- 
teracting fermions in one dimension, called a Tomonaga- 
Luttinger liquid, which is described by nonperturbative 
techniques such as bosonization [l|. For Landau- Fermi 
liquids the electron spectral function has a sharp peak 
with a nonzero weight, called the quasiparticle renormal- 
ization factor. The quasiparticle renormalization factor 
vanishes for interacting electrons in one dimension. 

For a FQHE state 0, the excitations in the bulk are 
gapped, but arbirtarily low-energy excitations are avail- 
able at the edge. The edge dynamics is one dimensional, 
and constitutes a realization of a chiral Luttinger liquid, 
which is a Luttinger liquid consisting of fermions moving 
only in one direction Here, only electrons at one 

edge of the FQHE system are considered; coupling with 
the oppositely moving electrons at the other edge is ne- 
glected, which is a good approximation for wide samples 
for which the two edges are spatially far separated. Wen 
has proposed an effective chiral Luttinger liquid (ECLL) 
model for the description of the FQHE edge , accord- 
ing to which the long distance properties at the edge are 
universal, described by a quantized exponent the value of 
which is determined by the quantized Hall resistance of 
the bulk FQHE state. The most direct probe of the prop- 
erties of this liquid is tunneling of an external electron 
laterally into the edge of the FQHE system. Non-linear 1- 
V characteristics have demonstrated NLFL behavior for 
the FQHE edge. The ECLL approach predicts an I-V 
behavior for the tunnel conductivity of the form / oc V"' , 
with a universal value for a. In particular, for FQHE 
states at filling factors v — n/{2pn + 1), the effective ap- 
proach predicts a ~ 2p+\. Experiments find a nontrivial 
value for a (i.e. a 7^ 1), but the observed a deviates from 



the predicted one HI 



. This discrepancy has motivated 



much work HiilUllil, including the present paper. 

The ECLL approach is built upon the idea of a one-to- 
one correspondence between the fermionic and bosonic 
Fock spaces in one dimension, and identifies a relation- 
ship between the operators of the two problems. In par- 
ticular, the fermionic field operator il}{x) is related to 
the bosonic field operator ^(x) through the expression 
ijj(x) ~ exp[— i(/()(a;)], which can be established rigorously 
for the ordinary one dimensional systems [l|. A similar 
rigorous derivation has not been possible for the elec- 
tron field operator at the edge of a FQHE system. Wen 
postulates that the electron operator at the edge of the 
I' = 1/m FQHE system is given by 



(1) 



which has the virtue of satisfying the antisymmetry prop- 
erty when m = 2p -|- 1 is an odd integer . This form 
leads to the quantized exponent for the I-V of the tun- 
nel conductance. It is not known at the present what 
causes the discrepancy between the effective theory and 
experiment. 

Our work is an outgrowth of the exact diagonalization 
studies of Palacios and MacDonald [l^ at = 1/3, the 
results of which were interpreted by the authors as con- 
firming Wen's ansatz. We now know, however, that the 
discrepancy is small at v = 1/3 (the observed value of 
a « 2.7 is close to the predicted a — i), and the results 
of Ref. [m were not sufficiently accurate to capture such 
small deviations. We extend their calculations to larger 
systems, using a Monte Carlo method, to obtain more 
accurate thermodynamic extrapolations. Furthermore, 
Ref. [1^ assumed a short range interaction model; we 
also investigate to what extent the results are sensitive 
to the form of the interaction. Our results are consistent 
with the ECLL predictions for the short-range model, 
but inconclusive for the Coulomb interaction. 

It should be stressed that the FQHE liquid itself 
provides a beautiful paradigm for a breakdown of the 
Landau-Fermi liquid concept. Here, strong interactions 
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generate electron-vortex bound states called composite 
fermions, which are qualitatively distinct from, and per- 
turbatively unrelated to, electrons Composite 
fermions experience a greatly reduced effective magnetic 
field, and possess quantum numbers (for the local charge 
and braiding statistics) which are a fraction of the elec- 
tron quantum numbers. The bulk properties of the com- 
posite fermion (CF) liquid have been investigated by a 
variety of means, and many experiments have directly 
verified the effective magnetic field, a clear indication 
of the NLFL nature of the state. The question of how 
an external electron, tunneling vertically in the CF liq- 
uid, couples to the system has been studied, and it has 
been predicted that it tunnels resonantly into an excited 
state, which is a bound state of several excited composite 
fermions, to produce a sharp peak in the electron spec- 
tral function 16lll7j . The CF liquid therefore provides 



a different mechanism for the breakdown of the Landau- 
Fermi liquid concept: The electron renormalization fac- 
tor remains nonzero, but lower energy states appear that 
are not described in terms of new quasiparticles. 

The plan of the paper is as follows. Section II lists 
the connection between the bosonic approach and wave 
functions. The results are given in Sec. Ill, followed by 
conclusion in Sec. IV. 



II. SPECTRAL WEIGHTS 

Following Palacios and MacDonald [l^ we test the va- 
lidity of Eq. [T] by comparing the microscopically calcu- 
lated spectral weights 



a 



(2) 



(where 9 denotes the position in one dimension, wrapped 
into a circle via periodic boundary conditions) with the 
predictions of the ECLL approach. In the ECLL ap- 
proach at filling factor 1/m the vacuum state |0) 
contains no bosons and the various symbols have the fol- 
lowing meaning: 



' ze 



(3) 



(4) 



(5) 



and 



\{ni})=Y[ 



1=0 



'nil 



-|0). 



(6) 



Here a| and a/ are creation and annihilation operators 
for a boson in the angular momentum I state, with the 
total angular momentum given by 



AM 



J2lni. 



(7) 



By expanding exp[i-ym(/)+(6')] and exp[—iy/rn(j)-{d)] in 
power series of products of a/ and a|, it is straightfor- 
ward to compute the analytical values for the spectral 
weights C{„j. The square of the spectral weight |C{„,}p 
is independent of the angular position parameter. The 
denominator in Eq. [2] eliminates the unknown normal- 
ization constant in Eq. [3] 

From the perspective of electrons, the vacuum state is 
the ground state of interacting electrons at v = l/m, and 
the field operator has the standard meaning of 



(8) 



where cj and ci are creation and annihilation operators 
for an electron in the angular momentum I state, the 
wave function for which is rji : 



m{z) 



z = X — ly. 



V2tt2H\ 

The denominator of Eq. [2] is interpreted as 

(01^^^0)10) = «+i|V'^(0)lO 



(9) 



(10) 



where VP^ is the normalized ground state of N interacting 
electrons at = l/m. The numerator is interpreted as 

({ni}\4>\em = (*f4i^l(f?)|0. (11) 



Here, we have, 

(12) 

where A is the antisymmetrization operator and J\fc is 
the normalization constant. 



The 



c function ^'i , , the electronic counterpart 
of the bosonic state represents an excited state 

that involves increasing the angular momentum of the 
N particle ground state by L units. It is not immedi- 
ately obvious how to construct it for a general case. For 
the integral quantum Hall state at v = 1, Stone showed 
[IMiillo] that it is obtained by multiplying the (TV + 1 
particle) ground state by the factor Y[i S";"' , where Si are 
defined as: 



N 



(13) 



The product increases the total angular momentum of 
the A'' -I- 1 particle ground state by AM = '}2,ilni. 
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The composite-fermion analogy suggests the identifica- 
tion [H 

i*f„t})=-^»' (^n(^'r)< (14) 

at V — 1/m, where A/'„, is the normaUzation constant. 
The angular momentum changes L and AM, defined rel- 
ative to the N and TV -|- 1 particle ground states, respec- 
tively, are related by 

L = AM-fm7V, (15) 

where the last term is the difi^erence between the angu- 
lar momenta of the ground states. We will assume the 
identification in Eq. [T3]in what follows; further justifica- 
tion for it is given in the Appendix. In our calculations 
with exact diagonalization method, we will use that the 
second quantization representation for the operators Si 
is given by (apart from a constant factor) 

^'-f:#i^4..c.. (16) 

fe=0 

III. RESULTS 
A. Exact diagonalization: Short range interaction 

Palacios and MacDonald [l^ computed the squared 
spectral weights |C{„;}p defined in Eq. [2] by exact diag- 
onalization for a short-range interaction model for which 
the Laughlin wave function [2T'| and the excited states in 
Eq. [M] are exact eigenstates. (This interaction takes a 
nonzero value for the Vi pseudopotential [l^ but sets all 
other pseudopotentials to zero.) They obtain results for 
systems containing up to eight electrons for AA/ =1 to 4; 
from an extrapolation to the thermodynamic limit, they 
find approximate consistency with the predictions of the 
effective ECLL approach. 

The results we obtained for Vi interaction are shown in 
Table HI and are identical to the results by Palacios and 



MacDonald fl3^. However, our extrapolation to thermo- 
dynamic limit differs slightly from theirs. The extrapo- 
lation assumes a leading finite size correction for |C{„,}p 
proportional to 1/iV; this dependence has not been de- 
rived analytically, but matches the numerical data well 
for these particle numbers. The thermodynamic values 
are close to, but significantly different from, those pre- 
dicted by the ECLL model. 



B. "Exact" diagonalization: Coulomb interaction 

A proper extension of the above results to the Coulomb 
interaction is not known. In the above, both the ground 
and excited states are exact eigenstates of the Vi model. 
While the exact ground state for the Coulomb interaction 
can be obtained for small systems, no operators analo- 
gous to the Si of Eq. [T3]are known that produce exact ex- 
cited states. We use an approximate, "hybrid" approach, 
in which we take the exact Coulomb ground state, but 
use the same operators Si to create excited states. The 
above calculation can then be extended to the Coulomb 
interaction. The exact spectral weights with for angular 
momenta AM = 1 to 4 and particles N =4 to 9 are tabu- 
lated in Table HIl Again, there is a small, but significant, 
deviation from the ECLL results. 

We mention here a technical point relating to the ef- 
ficiency of the numerical calculation. The number of 
basis vectors increases rapidly with the number of par- 
ticles. For example, for AM = 0, TV = 8 and 9 and 
10, the number of basis vectors are 55,974, 403,016 and 
2,977,866, respectively. During bosonic state creation 
and overlap calculation, the most time consuming step 
is searching for a given component from the ket vector 
^'Il^cT) ^^'^ matching it with the corresponding bra vec- 
tor component (4';^^^ |. We collect the fermionic states 
|mi, 7712, . . . , ttin) into bins indexed by the smallest three 
angular momentum values {toi, TO2, ms}, which allows 
the search for a match to be restricted to a single bin. 
This reduces the computation time by a factor of 70, en- 
abling computation of single boson states within 35 hours 
on our supercomputing cluster. 



C. Monte Carlo simulation results 



The exact diagonalization method has the drawback 
of being restricted to small numbers of particles. Fortu- 
nately with quantum Monte Carlo methods we can ex- 
tend the results to much larger systems for the above- 
mentioned short-ranged Vi interaction, for which all the 
wave functions in question are explicitly known. With 



the Monte Carlo calculation method we have obtained 
the spectral weights for up to 40 particles, for AM rang- 
ing from 1 to 6. The results are shown in Table IIIII and 
Figs, fflandll 
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AM {m} TV = 4 iV = 5 TV = 6 TV = 7 TV = 8 oo ECLL 



{1000} 
{2000} 
{0100} 
{3000} 
{1100} 
{0010} 
{4000} 
{2100} 
{1010} 
{0200} 
{0001} 



2.6000 
3.6400 
1.2585 
3.6400 
3.6740 
0.9356 
2.9120 
5.6503 
2.7489 
1.0184 
0.9044 



2.6667 
3.7778 
1.2953 
3.7778 
3.8133 
0.9358 
2.9907 
5.8730 
2.7952 
1.0353 
0.8583 



2.7142 
3.8755 
1.3224 
3.8775 
3.9128 
0.9390 
3.0466 
6.0247 
2.8284 
1.0484 
0.8302 



2.7500 
3.9531 
1.3425 
3.9531 
3.9860 
0.9425 
3.088 
6.131 
2.852 
1.058 
0.8109 



2.778 


2.9532 


3 


4.012 


4.3769 


9/2 


1.358 


1.4562 


3/2 


4.012 


4.3777 


9/2 


4.041 


4.4053 


9/2 


0.946 


0.9543 


1 


3.121 


3.3267 


27/8 


6.209 


6.7701 


27/4 


2.869 


2.9894 


3 


1.064 


1.1102 


9/8 


0.797 


0.6881 


3/4 



TABLE I: Squared spectral weights |C{„j}|^ from exact diagonalization method for the short range model in which only the 
Vi pseudopotential is nonzero. is the number of electrons in the study. For finite N, our numbers are identical to those in 
Ref. [13I , but our extrapolations to A — > 00 are slightly difi'erent. "(A linear extrapolation in 1/A is assumed for obtaining the 
thermodynamic limits in this table, in Tables II and III, and in Figs. 1 and 2 for the exact diagonalization results; the validity 
of this assumption is questioned by larger system Monte Carlo results shown in Fig. 1 for the short-range interaction.) 



AM {m} A = 4A = 5A = 6A = 7A = 8A = 9N- 



ECLL 



{1000} 

2 {2000} 
{0100} 

3 {3000} 
{1100} 
{0010} 

4 {4000} 
{2100} 
{1010} 
{0200} 
{0001} 



2.6000 
3.6400 
1.2871 
3.6400 
3.7551 
1.0029 
2.9120 
5.7699 
2.9313 
1.0632 
0.9676 



2.6670 
3.7783 
1.3177 
3.7783 
3.8768 
0.9836 
2.9911 
5.9657 
2.9315 
1.0700 
0.9128 



2.7150 
3.8786 
1.3504 
3.8786 
3.9937 
0.9925 
3.0475 
6.1450 
2.9843 
1.0926 
0.9095 



2.7514 
3.9553 
1.3654 
3.9554 
4.0522 
0.9830 
3.0903 
6.2311 
2.9710 
1.0935 
0.8725 



2.7801 
4.0157 
1.3744 
4.0157 
4.0885 
0.9780 
3.1233 
6.2804 
2.9636 
1.0886 
0.8420 



2.8031 
4.0644 
1.3849 
4.0644 
4.1265 
0.9773 
3.1499 
6.3315 
2.9662 
1.0898 
0.8264 



2.9619 
4.3954 
1.4658 
4.3954 
4.4334 
0.9588 
3.3365 
6.8044 
3.0087 
1.1193 
0.7259 



3 
9/2 
3/2 
9/2 
9/2 

1 

27/8 
27/4 

3 
9/8 
3/4 



TABLE II: Squared spectral weights |C{„,}p obtained from exact diagonalization method for the Coulomb interaction. 



D. Orthogonality 



The ECLL approach also predicts orthogonality be- 
tween states |5'{„,}) and for {ni} ^ {n'l}- Such 
an orthogonality is not apparent from the wave function, 
and does not follow from any symmetry. We have numer- 
ically tested it for several cases; the results, summarized 
in Tabic llVi indicate that the overlap between different 
states rapidly diminishes with increasing number of par- 
ticles. This also demonstrates that the states generated 
by the operators Si do not exactly represent free boson 
states for finite N, but are meaningful only in the ther- 
modynamic limit. 



IV. DISCUSSION AND CONCLUSIONS 

The thermodynamic values of the spectral weights, 
obtained by assuming a linear fit as a function of 
1/A^, lie within a few percent of the ECLL pre- 
dictions (0.9 % to 18 % for single boson states; 
see Table but the deviations between the two 

are often significant. However, deviations from 
linear fit are seen for large N. For example, for 



{1000}, {0100}, {3000}, {1100}, {0100}, {0001}, {010100}, 
{00001}, {200100}, {1000010}, {01010}, {002000}, 
{000001} (Figs. □ and ID), the actual fit is nonlin- 
ear, and tends toward the ECLL results with increasing 
N. A proper thermodynamic value is difficult to 
estimate accurately in many cases, partly because of 
the lack of analytic results regarding the N dependence 
of |C{„j}p. However, the thermodynamic limits are in 
general agreement with the ECLL predictions. This is 
consistent with several previous studies 

that 

confirm the validity of the ECLL model for the short 
range interaction. 

With a partial inclusion of the Coulomb interaction, 
linear extrapolation of the results for up to nine particles 
again gives matrix elements that are close to the ECLL 
prediction but significantly different. The results are also 
in general quite different from those for short range in- 
teraction. Here, however, it is not possible to extend our 
study to larger systems. From our experience with the 
short range interaction, we cannot rule out deviation of 
the |C{„j}p vs. 1/A^ plot from linearity for large N, and 
therefore consider our study as being inconclusive. 

Our calculations essentially serve as a test of the iden- 
tification of the operators Si with the boson operators a/ 
for the short range interaction. This was first introduced 
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by Palacios and MacDonald, and is further justified in 
the appendix, but not proven rigorously. 
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AM {m} 


l^ini > 1 




% deviation 


1 


{1000} 


2 Q782('0 02^ 


3 


0.89 


2 


{2000} 


4 5136f0 02^ 


4.5 


0.30 




{0100} 


1 47finf0 01 


1.5 


1.60 


3 


{3000} 


4 4896f0 02) 


4.5 


0.23 




{1100} 


4 4764f0 02) 


4.5 


0.52 




{0010} 


9769(0 01) 


1 


2.31 


4 


{4000} 


3 3833^0 03) 


3.375 


0.25 




{2100} 


6 7928 fO 05) 


6.75 


0.63 




{1010} 


3 0005 fO 02) 


3 


0.02 




{0200} 


1.1273(0.01) 


1.125 


0.20 




{0001} 


0.7210(0.01) 


0.75 


3.87 


5 


{50000} 


2.0429(0.02) 


2.025 


0.88 




{31000} 


6.9042(0.04) 


6.75 


2.28 




{20100} 


4.5869(0.03) 


4.5 


1.93 




{10010} 


2.2493(0.02) 


2.25 


0.03 




{01100} 


1.5073(0.01) 


1.5 


0.49 




{00001} 


0.5481(0.01) 


0.6 


8.65 


6 


{600000} 


1.0303(0.01) 


1.0125 


1.76 




{410000} 


5.2688(0.04) 


5.0625 


4.08 




{301000} 


4.6980(0.04) 


4.5 


4.40 




{200100} 


3.5203(0.03) 


3.375 


4.31 




{220000} 


5.2251(0.04) 


5.0625 


3.21 




{100010} 


1.7837(0.02) 


1.8 


0.91 




{111000} 


4.7573(0.03) 


4.5 


5.72 




{010100} 


1.1669(0.01) 


1.125 


3.72 




{002000} 


0.5014(0.01) 


0.5 


0.28 




{000001} 


0.4072(0.01) 


0.5 


18.56 



TABLE IIL Monte Carlo simulation results for Vi pseudopotential with 20 particles. Third column shows the thermodynamic 
extrapolation of squared spectral weight. The statistical errors are shown in parenthesis. ECLL predictions and deviation from 
ECLL predictions are shown in fourth and fifth columns, respectively. 
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AM {m} {n'l} N^6 N^S N = 15 

~ {2000} {0100} 6.90(0.1) xlO"^ 4.3(0.12) xlO"^ 1.3(0.06) xlO"^ 

3 {3000} {1100} 2.1(0.05) xlO'^ 1.2(0.03) xlO"^ 4.0(0.06) xlO^^ 
{3000} {0010} 1.0(0.07) xlO""' 8.6(0.52) xlO-"' 1.2(0.57) xlO^^ 
{1100} {0010} 4.3(0.03) xlO-2 2.7(0.10) xlO"^ 7.7(0.24) xlO"^ 

4 {4000} {2100} 4.0(0.05) xlO'^ 2.5(0.03) xlO'^ 8.1(0.20) xlO"^ 
{4000} {1010} 4.20(0.2) xlO"-* 1.5(0.20) xlO"'' 4.0(1.00) xlO"^ 
{4000} {0200} 1.0(0.05) xlO^"' 5.0(0.30) xlO""^ 5.5(1.40) xlO"^ 
{4000} {0001} 4.20(1.2) xlO"'' 1.0(0.2.4) xlO"** 2.7(1.20) xlO"^ 
{2100} {1010} 8.1(0.05) xlO-2 5.0(0.04) xlO'^ 1.6(0.01) xlO'^ 
{2100} {0200} 1.4(0.02) xlO'^ 8.3(0.80) xlO"^ 2.7(0.10) xlO"^ 
{2100} {0001} 1.8(0.05) xlO-3 6.8(0.30) xlO"" 8.0(1.10) xlO'^ 
{1010} {0200} 1.1(0.03) xlO-^ 4.0(0.10) xlO"'' 4.2(0.90) xlO"^ 
{1010} {0001} 8.8(0.05) xlO-2 5.3(0.04) xlO'^ 1.6(0.02) xlO'^ 
{0200} {0001} 6.0(0.04) xlO~^ 3.5(0.03) xlO"^ 1.0(0.03) xlO"^ 



TABLE IV: Overlaps l*^^/})]^, where *{nj} is defined as in Eq. (fT4)) . Results are shown for AM=1 to 4 for several 

values of A''. 
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FIG. 1: (Color online) Squared spectral weights |C{„,}p obtained from several methods: exact diagonalization for short 
range interaction (labeled Vl-ExD) and Coulomb interaction (Coulomb-ExD), and Monte Carlo evaluation for the short range 
interaction (Vl-MC). The symbols are defined on the figures. The error bars indicate the statistical uncertainty in the Monte 
Carlo results. Results are shown for AM = 1, 2, 3, 4 and 5. For AM — 5, {50000} state we have only Monte Carlo results. The 
solid line and the dashed line are a linear fit for the exact diagonalization results for short range and Coulomb interactions, 
respectively. The chiral Luttinger liquid (CLL) predictions are marked by solid diamonds on the vertical axes. 
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FIG. 2: (Color online) Filled squares represent squared spectral weights |C{„,}p obtained from Monte Carlo simulaton for Vi 
interaction for AM = 5, 6. The error bars indicate the statistical uncertainty in the Monte Carlo results. Thermodynamic 
extrapolation is also shown as a blue square on the vertical axis. Their numerical values are given in TableHni ECLL predictions 
are shown as solid diamonds on the vertical axes. 
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APPENDIX A: BOSONIC OPERATORS 

In the standard bosonization approach for one- 
dimensional Fermi systems, the density operators play 
the role of bosons. For the state at i/ = 1, Stone ob- 
served that the fermionic excitations at the edge can 
be mapped into linear combinations of products of sym- 
metric polynomials Si 's given in Eq. (|13l ) : 



Ai+Af-l X2+N~2 

^1 ^1 

Ai+JV-1 \2+N-2 

^2 ^2 



-1 

-,AiY 



Ai+AT-l X2+N-2 Ajv 



N 

{0}' 



(Al) 



where is the ground state at ly 



1, C}^\^'s are the 



expansion coefhcients (cf. Eq 4.12 in [18|), Ai > A2 > 
• • • Aat, and J^k '^fc ~ ^'^i = is the total angular 
momentum of the excited state, measured relative to the 
ground state. The symmetric polynomial Si is identified 
with bosonic operator aj at angular momentum I. We 
now ask under what conditions the operators Si satisfy 
the canonical commutation relations. For this purpose, 
we use the form given in Eq. (jl6p . 

The commutation relation [Sk,Si] — 0, where k and 
I denote angular momentum quantum numbers, follows 
from the fact that the S"s can be expressed as polynomi- 
als, and can also be verified straightforwardly using the 
operator form of S in Eq. (|16p . Demonstration of 



(A2) 



[Sk,Sj] cx Ski 
is more tricky. An explicit evaluation gives 



{m + ky.{m + ly. I 



ml(m + k-iy. 



the right hand side of which, in general, does not vanish 
for k^l. 

To proceed further, we assume the idealized occupation 
number for the ground state 



{n\clck\n) = 



if fc < M 
otherwise , 



where M is the angular momentum of the outermost oc- 
cupied orbital. This is surely correct for k << M, but 
only approximate near the edge. We further assume that 
the effect of the operators Si is confined to the edge; this 
is manifestly correct for v = I, but not obvious for the 
FQHE state = 1/m. 

We first consider k = I. Defining normal ordering of 
operators in the usual manner (i.e., by subtracting the 



ground state expectation value), we get 
(to -|- fc)! 



[<5'fe, S"; 



E 



Cm Cm 



+ {n\cl,^ck+r,M-{^\clck\n) 



(A4) 



The dominant contribution comes from the last two 
terms on the right hand side. (When only terms near 
the edge contribute to the commutator, then for to « M, 
the factorial term varies as M*^, and the normal ordered 
terms cancel to the lowest order.) This gives 



(m + fc)! 



m=0 



TO! 



[{n\cl^^ck+^\n) - {n\clck\n) 



M 



E 

m=M-fe+l 

-vkM''. 



i^^(0-.) 
to! 



(A5) 



The canonical bosonic relations are obtained by defining 

1 



4 = 



ViJkM^ 



Sk, fc>o. 



(A6) 



Next we consider the case when k ^ I, assuming fc > I 
and N » {fc, 1} (since only excitations near the edge are 
significant). The commutator can be expressed as 



[Sk,Sj] 



E 



' (m + fc — ly.ml 
- 7V2 



(to — I) 



I (m + fc)!2 
to,! (to, + fc — ^)! 



^'rn+k-l^- 



(A7) 



We have used above that the vacuum expectation values 
of all combinations on the right hand side vanish due to 
angular momentum conservation. Only terms for which 
TO « M are significant in the summation. A little algebra 
using the asymptotic expansion 



T{n p) 
r(n + q) 



1 



{p-q){,p + q- 1) 



12n2 



2n 

p-q 



{i{p + q-lf -p + q-l) 
(A8) 



shows that the commutator vanishes to leading order, 
producing 



[Sk,S\]^ ^ V^^^O(1/to2)c^ 



m+k-l'^rn- 



(A9) 



For the bosonic operators defined in Eq. (|A6p . this im- 
plies 



E V(^)'^'^(l/m')cL+fe_,c™-.0, 

(AlO) 
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thus completing the relationship between Si and bj . We but assume that the effect of the operators Si is confined 
note that these considerations are valid for arbitrary i^, to the edge. 
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